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We present a comprehensive experimental and theoretical studies of photoluminescence of single 
Cdi_2,Mn2,Te quantum dots with Mn content x ranging from 0.01 to 0.2. We distinguish three 
stages of the equilibration of the exciton-Mn ion spin system and show that the intermediate stage, 
in which the exciton spin is relaxed, while the total equilibrium is not attained, gives rise to a specific 
asymmetric shape of the photoluminescence spectrum. From an excellent agreement between the 
measured and calculated spectra we are able to evaluate the exciton localization volume, number 
of paramagnetic Mn ions, and their temperature for each particular dot. We discuss the values of 
these parameters and compare them with results of other experiments. Furthermore, we analyze 
the dependence of average Zeeman shifts and transition linewidths on the Mn content and point out 
specific processes, which control these values at particular Mn concentrations. 

PACS numbers: 78.67.Hc 78.55.Et 71.70.Gm 75.75.-c 
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I. INTRODUCTION 

Zero dimensional density of states for carriers con- 
fined in self-assembled quantum dots (QDs), make these 
nanostructures unique in many aspects. Most impor- 
tantly from the point of view of applications, they of- 
fer a combination of the zero-dimensional, atomic-like 
optical properties with easy integration into the exist- 
ing semiconductor technology.^ As a result, many pro- 
posals emerge on employing quantum dots as sources of 
classical and non-classical light, El quantum bits (qubits) 
and quantum gates,** ^ systems for storage of classical^ 
and quantunt^ information, and ultra-precise sensors of 
local electric fields .El 

Semimagnetic QDs possess an additional degree of 
freedom i.e., the magnetic moment associated with the 
magnetic impurities incorporated into the dot material. 
Therefore, these dots offer a unique opportunity to com- 
bine the zero-dimensional electronic properties with the 
paramagnetism of the semimagnetic semiconductor. Har- 
nessing the magnetism in these nanostructures can lead 
to novel nanoscale devices. In particular, electrical con- 
trol of the QD charge state is expected to result in switch- 
ing of the QD magnetization on and off.ll^^ Optical con- 
trol of the magnetizatiorPS and its long relaxation time^ 
can provide a path toward employing semimagnetic dots 
in memory storage devices. Also, a pair of coupled non- 
magnetic and semimagnetic dots was proposed as a qubit 
gate, with electric field tunability.ES 

By far, the most studied semimagnetic QDs are made 
of II - VI semiconductors, since the low solubility of Mn 
in III - V compounds hinders introduction of more than 
one Mn ion into a dot.*® Among the studied materials, 
the most attention has been devoted to Cdi_j;Mnj;Te and 
Cdi_a;Mna;Se. In these dots, the s,p-d exchange interac- 
tion between the charge carriers and localized magnetic 



ions underlies the optical properties analogously to higher 
dimension systems .1^^ In particular, giant Zeeman split- 
tings of the excitonic sublevels are observed, with a mag- 
nitude that is directly proportional to the dot magnetiza- 
tion as evidenced by quantitative agreement between the 
exciton recombination energies and a modified Brillouin 
function.liSHSfll Another striking spectroscopic feature of 
these dots is a substantial broa dening of the photolu- 
minescence (PL) transition .ESiSl xhe linewidths at zero 
magnetic field are on the order of a few meV, whereas for 
nonmagnetic CdTe and InAs QDs, the transitions ar e one 
and two orders of magnitude narrower, respectively P* ! ^^ ! 
The broadening is another consequence of the exciton- 
Mn exchange interaction: thermal fiuctuations of the Mn 
spins introduce a temporal magnetic field, fiuctuating on 
the timescale of hundreds of nanoseconds, which in turn 
infiuences the transition linewidth in a time averaged 
spectrum. Remarkably, with increasing magnetic field, 
the PL transitions become more narrow, as the fiuctua- 
tions of the Mn spins are suppressed when their Zeeman 
splitting becomes larger then the thermal energy. Thus, 
at high magnetic fields, the transition linewidth tends 
to the one observed for a Mn-undoped, CdTe or CdSe 

dot.cmnmi 

Most of the studies on semimagnetic dots involved 
an assumption of thermal equilibrium between the Mn 
spin system and the exciton. In such a case, the ex- 
change energy is minimized by spontaneous formation 
of a fe rromagn etic order in the paramagnetic Mn spin 
system.liSiSIMMI This spontaneous magnetization forms 
within the exciton wave function^^ and is hence labeled 
as excitonic magnetic polaron (EMP). Formation of the 
EMP is therefore synonymous with reaching the equi- 
librium in the exciton-Mn spin system. The lineshapes 
become Gaussian with the full width at half maximum 
(FWHM) directly monitoring the statistical spin fiuc- 



tuations in the exchange field of the excitonllSl How- 
ever, achieving the equilibrium requires the character- 
istic tiniescale of collective Mn magnetization dynam- 
ics, closely related to Mn spin-spin relaxation time,'^*' to 
be shorter than the exciton recombination time r^. In 
Cdi_a;Mna;Te QDs, the latter time is on the order of 300 
psP21 Thus, in order to achieve the equilibrium, the Mn 
magnetization dynamics has to occur on a tiniescale of 
1-100 ps, much shorter than those driven by spin-l attice 
interactions, which occur on the /xs/ms timescales P" * ^^ ! 
The spin-spin mechanism becomes efficient only if the av- 
erage distance between Mn spins is small enough. Con- 
sequently, EMP formation was unambiguously demon- 
strated only for dots with a Mn content x > 0.1. For 
dots with a more diluted concentration of Mn spins, EMP 
starts to form upon photoexcitation, but exciton recom- 
bination interrupts the establishment of equilibrium.'^^l 

Another important question con cerns the possibility of 
an optical orientation of the EMP.^iSEll This could only 
be possible if the exciton spin relaxation time were longer 
than the EMP formation time Tf. We have shown that 
for relatively large self-assembled Cdi_a:Mn2:Te QDs the 
opposite is true.^^ We have established the EMP forma- 
tion scenario in which the exciton adjusts its spin to the 
direction of a momentary, initial magnetization fluctu- 
ation in the dot and then, provided that r^ > Tf the 
EMP is formed. Therefore, we distinguish three stages 
in reaching the equilibrium. Initially, shortly after the 
photocreation of the exciton, the system is completely 
unrelaxed with the exciton spin being either random if 
the orientation was lost during the capturing by the QD 
potential, or reflecting the photon angular momentum if 
otherwise. In either case, the exciton spin has no prefer- 
ential orientation with respect to the magnetization. In 
the second stage, the exciton relaxes its spin, which is 
synonymous with thermal population of the exciton spin 
states in the exchange field of the Mn spins. The Mn 
magnetization is then increasing due to Mn spin-spin in- 
teractions until it reaches its equilibrium value for a given 
exchange field provided by a spin-polarized exciton. At 
the end of this third stage the full equilibrium is reached 
and the EMP is formed. 

In this paper, we aim to establish a quantitative de- 
scription of the optical properties of semimagnetic QDs 
containing many Mn ions and demonstrate the influence 
of the three relaxation stages on the shape of the PL 
spectrum. In particular, we show that the relaxation of 
the exciton spin has a nontrivial effect on the shape of 
the PL line and leads to its narrowing by a factor as 
high as 2 and development of asymmetrical shape in a 
certain range of temperatures and magnetic fields. We 
start with presenting the theoretical model, which allows 
to compute the PL spectra particular for the different re- 
laxation stages. Among model parameters are the num- 
ber of paramagnetic Mn ions interacting with the exciton 
and the exciton localization volume, related to the QD 
volume. Thus, the by fltting the model to the measured 
PL spectra, we gain access to vital structural parameters 



describing these semimagnetic dots. We start present- 
ing the experimental results with showing how the PL 
spectrum evolves in time reaching subsequent stages and 
confirming our predictions. We then demonstrate the 
change of the line shape as the magnetic field and tem- 
perature is varied. We directly compare the measured 
spectra with the calculated ones and find an excellent 
agreement, which allows us to estimate the structural 
parameters mentioned above. Furthermore, we analyze 
Zeeman shifts and PL line widths for a large set of QDs 
from samples with the Mn content varied between 0.01 
and 0.2 and discuss them in the framework of the pre- 
sented model. 



II. THEORETICAL DESCRIPTION 

A. Distribution of Mn spins in an external field 

In this work, we are interested in statistical properties 
of a finite number of Mn spins (A^Mn) confined in the QD 
volume. For a given Mn concentration x, A^Mn experi- 
ences statistical fluctuations with a standard deviation 
SNmh oc ^/Nmu, making it more meaningful to use A^Mn 
instead of x as an independent variable. We study QDs 
with A^Mn > 20 and at moderate magnetic fields, where 
the magnetization is far from saturation. This allows us 
to treat the total Mn spin S = J^i ^'" S^ as a classical 
variable and assume that the probability distribution of 
S^ is Gaussian: 



p{sn = 



1 



v^ 



: exp 



(S,^,. „ 



'S \ "^s 



In the above, the variance of S^, ag is given by: 



a| = ((^-)2) - (5^)2 = -kgT 
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gfiB OB 



{S' 
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which follows from application of the Fluctuation- 
Dissipation Theorem with T being the Mn spin temper- 
ature. The magnetic field dependence of (S^) is given by 
the well-established, phenomenologically modified Bril- 
louin function,l32l 



(5^ 



-A^Mn ■ S ■ Bs 



gfj'B ■ s 

ksiT + To) 
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where g, fiSi and ks are the Mn ion g— factor, the 
Bohr magneton, and Boltzmann constant, respectively, 
and the spin of an individual Mn ion S — 5/2. The 
Mn impurity spins in a II- VI material can be considered 
paramagnetic only in the limit of very small concentra- 
tion X. At X > 0.01 the short-range antiferromagnetic 
interactions (both isotropic and Dzyaloshynsky-Moriya 
anisotropic ones) between the Mn spins have to be taken 
into account. Since in this work we are interested in QDs 
with a rather large Mn contents, A^Mn is usually sub- 
stantially smaller than the number of Mn ions physically 



present within the QD volume V. This is because the 
Mn spins located on nearest neighbor cation sites are an- 
tiferroniagnetically coupled with isotropic superexchange 
(with interaction strength J « 1 meV), thus removing 
these spins from paramagnetic response for fields smaller 
than « 10 T. In the following discussion of optical proper- 
ties of the Cdi_2.Mn2,Te QDs, we will only consider these 
Mn ions, which do not possess a Mn as a nearest neigh- 
bor and label them as paramagnetic Mn ions, which is 
not strictly correct since next nearest neighbors are also 
antiferromagnetically coupled and also decrease the para- 
magnetic response to the magnetic field. Therefore, the 
number of paramagnetic Mn ions Nmu will be treated 
as a model parameter accounting for all the Mn ions re- 
sponding to fields below 10 T. In the expression (Is]), the 
replacement of temperature T with T + Tq (x) accounts 
for a suppression of paramagnetic response due to the 
farther-neighbor Mn-Mn couplings. In this way we in- 
corporate the realistic description of interacting Mn (i.e., 
only approximately paramagnetic) spins into the calcu- 
lation of various statistical properties of the exciton-Mn 
complex.EalMll 



B. Description of exchange-coupled system of Mn 
spins and an exciton 

The Hamiltonian of the electron-hole pair (the exciton) 
and the Mn spins is given by: 



H ^ H, 



H, 



p—d 



H. + H 



X 



(4) 



consisting of the s-d (electron-Mn) exchange interaction, 
and the p-d (hole-Mn) exchange interaction, the Zeeman 
terms, and the Hx term related to the electron-hole inter- 
action within the exciton. The sp-d exchange interaction 
is given by^'' 

Hs-d + Hp^d = 

- a ^ <5(r, - ROS, • s - ^ ^ 6{vn - R.)S. • j/3 , (5) 

i i 

where a (/3) is the exchange constant for electrons (holes) 
with iVoa = 0.22 eV (A^o/3 = -0.88 eV) in Cdi.^jMn^Te, r^ 
(r/i) is the position operator of an electron (a hole), s (j) 
is the spin operator of an electron (a hole) , R^ is the posi- 
tion and Si is the spin operator of the i-th Mn ion. A^o is 
the cation density. We consider an exciton in its orbital 
ground state, i.e., both the electron and the hole occupy- 
ing respective lowest-energy orbitals, and we assume that 
its wavefunction is a product of electron and hole wave- 
functions. We note that this is only true for strongly con- 
fined excitons, where Coulomb correlation effects can be 
neglected. We obtain the Hamiltonian of the exciton-Mn 
exchange by projecting Hx-Mn = Hs-d+Hp-d onto two, 
two-dimensional subspaces corresponding to the electron 
and the hole. For the electron, we obtain the s-d exchange 



Hf,^n 



-a^|*e(R.)|'S, 



(6) 



where ^e(r) is the electron envelope wavefunction. Such 
a projection leads to a more complicated result in the 
case of holes. In low-dimensional structures such as QDs, 
the quantum confinement and strain splits the heavy hole 
(hh) and light hole (Ih) states. In CdTe QDs, the topmost 
valence band orbital is a twofold spin-degenerate state 
of mostly hh character, and the projection of the p-d 
exchange operator on this state givesPS 



^P-d = -/?El^''(^*)l'^^'j"/3- 



(7) 



where ^/i(r) is the hole envelope wavefunction. We 
note that keeping a finite admixture of Ih state in the 
topmost valence orbital changes the S^j^/3 term to 
{S^K^ +eS^K^ +e*5^K+)/2, where Ki are the Pauli op- 
erators defined in the two-dimensional subspace of hole 
orbital ground state, and |e| ^ 1 is a measure of the 
amount of Ih admixture.'^SHUl xhe mixing between heavy 
and light hole states in CdTe dots is much stronger than 
in InGaAs dots, resulting in such effects as linear polar- 
ization of the trion PL ,1221 non-zero in-plane holeo- factor 
and in-plane emission from the dark excitons,'^. How- 
ever, for the effects considered here it is of minor impor- 
tance, and in the calculations presented in this paper we 
will neglect it (i.e., we put e = 0). 

For magnetic field along the z— axis, the Zeeman term 
is given by: 



Hz = gefJ-BBs"" + gh^BBf + g^isS'' 



(8) 



where s^ is the z component of the electron spin operator, 
j^ is the z component of the hole angular momentum 
operator, and S^ corresponds to the total Mn spin, i.e.. 

The exciton Hamiltonian Hx includes both isotropic 
and anisotropic exchange interactions.^*^ The energy scale 
of the latter term is ~0.1 meV and thus we neglect it as 
it is some two orders of magnitude smaller than typical 
sp-d exchange energies considered here. The former term 
is: 



Hx = 



So 



11 



1-1) (-11 -|2) (21-1-2) (-21) 



(9) 
in which (Jq ~ 1 meV is the isotropic e-h exchange split- 
ting, |±1) are the bright exciton states (with electron 
and hole spins antiparallel) optically active in a± polar- 
ization, and |±2) are the dark exciton states (with elec- 
tron and hole spins parallel) . The recombination of dark 
excitons, as mentioned above, is allowed by the mixing 
between the hh and Ih subbands, but this mechanism 
is of negligible importance in Mn-doped QDs since here 
the main mechanism of bright-dar k exci ton mixing is due 
to fiip-fiop terms of sp-d exchange^ZEH. Below, we show 



that the influence of the dark states on the PL Hne shape 
of the dots considered here (i.e., containing at least 20 
Mn ions) is negligible. Thus, we will neglect Hx in the 
following calculations. The only indirect effect of Hx 
which we keep, is that together with Hp-d it suppresses 
the electron-Mn spin flip-flops (i.e. the S^s^ terms in 



Eq. ([6|)). It allows us to write the exciton-Mn exchange 
interaction neglecting the off-diagonal terms in S^ . Thus, 
we get: 



H 



X-Mn 



J2S!H'^e{R^ 



|2„z 



■/?|vI/,(R,)|V/3] 



(10) 
The final approximation that we will make is the as- 
sumption that the wavefunctions of both an electron and 
a hole are of a box (muffin-tin) shape, i.e., the elec- 
tron and hole wave functions are approximated by a con- 
stant value inside a volume Vg/^: ^e^/i(r) = 



e/h- 

while vanishing outside. This approximation, although 
crude, allows for easy diagonalization of Hx-Mn and was 
shown to correctly repro duce the static properties of the 
exciton-Mn system .l^SHI] In the muffin-tin approximation, 
the exciton-Mn interaction reads: 



H 



X-Mn 



= -s'- 



K 






^'73 



(11) 



C. Modeling of the PL spectra 

The two distinct confinement volumes for the electron 
and the hole in Hamiltonian (fill enable us to treat sep- 
arately the confinement of both carriers, a priori un- 
known. On the one hand, negligible band offset in the 
CdTe/ZnTe interface suggests a weakly bound hole, lead- 
ing to a much stronger leakage of the hole wavefunction 
out of the dot. On the other hand, PL experiments on 
CdTe dots with single Mn ions allow to estimate the rel- 
ative electron and hole confinement a nd po int out that 
the hole is in fact more strongly bound.l^ZISl xhe electron 
confinement volume is usually identified with the volume 
of the QD. In the following, we define an effective con- 
finement volume V . Taking into account the ratio of 
|/3/a| =4 and the realistic bound of 1/2 < Vh/Ve <2, one 
can show that for both cases: (i) Mn spins located only 
within the QD, and (ii) the same Mn concentration in 
the barrier and the dot, the effective confinement volume 
y is to a very good approximation equal to the volume 
occupied by the hole (see Appendix [A| for details). Using 
the Hamiltonian (111 for a given value of S^- we obtain 
the energies of the bright |±1) excitons. 

^^ ' ^ 2\Vh Ve) 2 2 V 2 ' 

(12) 
where Ez is the Zeeman splitting of the bright exciton 
states due to applied magnetic field. 

The PL intensity in a± polarization corresponding to 
recombination of |±1) excitons at energy E relative to 



the energy of a bare exciton (noninteracting with the Mn 
ions, its energy determined solely by the QD morphology) 
is then given by: 

I±{E) ^ j P{S')p{E±{S'))5{E - E±{S'))dS'- , (13) 

where p{E±) is the probability of occupation of exciton 
|±1) states at the moment of recombination. Note that 
we assume that the nonresonant excitation leads to cre- 
ation of both |±1) with equal probability. 

In the completely unrelaxed case we have p{E±) — 1/2 
and the exciton just probes the statistical distribution of 
the Mn ions. The PL spectrum is then given by: 



I±{E) ex exp 



[E T Eof 
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where as is given by Eq. (pi, and 
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(14) 



(15) 



(16) 



where (S^) is given by Eq. (pi. Note that in a nonzero 
B field the exciton Zeeman splitting Ez is typically neg- 
ligible compared to the first term in the above equation, 
and in the following we neglect the presence of Ez term. 
We see that in the completely unrelaxed case the line 
shape is Gaussian with the energy position and FWHM 
(7) determined by the mean and variance of Sz, respec- 
tively, which in turn are magnetic field dependent. Us- 
ing the above expressions FWHM of the PL transition is 
given by: 
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.... (1^) 
where the argument of the Brillouin function Bg is the 

same as in Eq. [s] In particular, at -B = 0, (Sz) — 
and the PL transition is centered around the energy of a 
bare exciton (see Fig. IT] yellow (light) curve). Naturally, 
the expression ( [T4| reproduces both the splitting and the 
narrowing of the transition line with increasing magnetic 
field (see line in Fig. [2[a)). 

If the exciton spin is relaxed, the occupation of the 
|±1) exciton states split by the sp-d exchange interac- 
tion is thermal and the majority spin population at low 
enough temperature probes roughly a half of the Mn spin 
states — those aligned antiparallel with the majority ex- 
citon spin. We take this into account by introducing the 
occupation factors: 



PiE±iSn) 



exp( 



ksTx ' 



2cosh( 



-E±(S- 
knTx 



(18) 



where we have again neglected the Ez term. In the 
above, we have introduced a temperature Tx, which cor- 
responds to the exciton spin temperature in the limit 
of the exciton spin relaxation time much shorter than 
recombination time. Since the exciton relaxation time 
decreases with increasing Mn concentration,'^^, this is 
true only for sufficiently large x and thus in principle 
Tx represents an upper limit for the exact exciton spin 
temperature. Indeed, in a tirsie- unresolved spectrum, we 
integrate both the recombination of spin-unrelaxed and 
spin-relaxed excitons. Quantitatively, Tx characterizes 
the degree to which the exciton spin population has equi- 
librated with the reservoir to which the exciton spin de- 
gree of freedom is most strongly coupled. Using Eq. ( 13 ) 
we get: 



I±{E) ex exp 



2a| 



2cosh(^) 



(19) 



An immediate consequence of this expression is that with 
increasing positive magnetic field (i? > leading to 
(5"^) < and Eq < 0) we obtain a redshift of cr+ emis- 
sion, while the cr_ emission is blueshifted, and its inten- 
sity drops rapidly to zero once liSo/fesTxl becomes larger 
than one. However, a more interesting consequence is 
that in the case of relaxed exciton spin, at zero magnetic 
field, the PL line shape is no longer Gaussian, but be- 
comes suppressed from the high energy side due to the 
Bolztmann factor (see Fig. [I] green (dark) curve). With 
decreasing the temperature, the asymmetry is enhanced 
as a result of larger imbalance in the thermal occupation 
of the exciton spin states (see Fig. l2jb)). On the other 
hand, as the magnetic field is increased, the exchange- 
driven spin splitting of the exciton states increases and 
the occupation probability of the spin-up state quickly 
approaches unity. Thus, the PL lines recover the Gaus- 
sian shapes for magnetic fields of about 1 T depending 
on temperature (see Fig. |2ja)). 

In order to model the PL spectrum for the exciton- 
Mn spin system in full thermal equilibrium, we have to 
take into account the formation of the polaron (see In- 
troduction) . The polarization of the Mn ions occurs as a 
result of an effective magnetic field B^x imposed by the 
exchange interaction with the exciton: 



B, 



2gtiBV 



(20) 



Thus, once the Mn population achieves a thermal equi- 
librium in this field, a net magnetization is created, 
and a redshift of the PL peak is developed. Impor- 
tantly, at external magnetic field Bq — the orienta- 
tion of the developed magnetization is random, as it de- 
pends on the random initial magnetization of fiuctuating 
Mn spin ensemble ,'23 resulting in no net polarization of 
PL from the redshifted peak (unless a memory effect is 
presenlP^. With increasing Bq, the PL rapidly becomes 
tT+-polarized. The polarization-unresolved PL signal in 
the regime of polaron formation is thus given by Eq. ( 19 ) 



albeit with Eq and cte calculated at the value of the total 
field of B — Bq + Bexf^ Since the effective field Bex is on 
the order of a few tesla leading to a redshift of the PL 
peak exceeding aE, we expect the lineshape for a fully 
equilibrated system to be a Gaussian. 
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FIG. 1. Calculated PL spectra for a QD with an exciton lo- 
cahzation volume covering NqV = 4000 cation sites and with 
A'^Mn = 40 paramagnetic Mn ions at a temperature of 4 K. 
Three stages of the relaxation process are modeled. Points 
demonstrate result of numerical diagonalization of a Haniil- 
tonian, which includes the off diagonal terms, neglected in 
developing the analytical model presented in the text. Lines 
denote the results of the analytical model. 
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FIG. 2. Magnetic field (a) and temperature (b) dependence 
of calculated transition linewidths (FWHM) for a QD with 
an exciton localization volume covering A^oV' = 4000 cation 
sites and with Ai^Mn = 40 Mn ions. Points (lines) denote the 
linewidths s of a recombination of fully spin-relaxed (com- 
pletely unrelaxed) excitons. 

In Fig. [T] we demonstrate the results of our model 
calculations for a dot confining an exciton to a volume 
spanning over NqV = 4000 cation sites with A^Mn = 40 
Mn ions at a temperature of 4 K. Points denote the 
PL spectrum obtained numerically by diagonalization of 
the Hamilt onian, which includes the off-diagonal spin-ffip 
termd^MU ^nd assumes full thermal equilibrium. This 
result is compared to the spectrum (dark blue line) cal- 
culated with an analytical formula (14 1 at B^x = 1-2 



T. which corresponds to the assumed locaHzation vol- 
ume V (see Eq. (20 1). We note that the agreement, 



despite dropping the off-diagonal terms in developing ex- 
pression (14 1, is very good. Also, for a fully thermalized 



exciton-Mn ions system, the spectrum is symmetric and 
Gaussian, as expected. The redshift of the exciton energy 
directly reflects the gain of the energy of the system upon 
the build-up of the magnetization. In this case, the mean 
exciton energy is more than 5 meV lower than for a dot 
where the full equilibrium was not reached. The spec- 
trum calculated for the unrelaxed case is just a Gaussian 
centered at zero — the energy of a bare, noninteracting 
exciton. However, for the spin-relaxed case, we find that 
the spectrum is strongly asymmetric and the linewidth is 
substantially reduced with respect to both the unrelaxed 
and the fully relaxed case. This is a strong indication that 
in the case of spin-relaxed excitons, the linewidth depen- 
dence on magnetic field can be non-monotonic. Indeed, 
as shown in Fig. l2Fa), the PL transition linewidth consid- 
ered in Fig. [TJincreases with the magnetic field up to 1 T 
and then decreases and follows the dependence for unre- 
laxed excitons, i.e., directly monitors the magnetization 
fluctuations, which are then suppressed with increasing 
the magnetic field. In Fig. p[b), we compare the temper- 
ature dependence of the linewidth for the unrelaxed and 
spin-relaxed cases. It demonstrates that the narrowing 
due to exciton spin relaxation is temperature dependent 
and that at liquid helium temperatures the linewidth can 
be twice smaller than in the unrelaxed case. 



III. SAMPLES AND EXPERIMENT 

The samples are grown by molecular beam epitaxy 
on (lOO)-oriented GaAs substrates. First, a 4 /im-thick 
CdTe buffer layer is deposited. Then, the first barrier 
layer, about 0.6 /um thick is grown. Cdi_a;Mna;Te dots 
are formed on top of this layer and covered with another 
barrier layer, 50 nm thick. QDs are formed using the tel- 
lurium desorption methocP^Bl from a 6 monolayer-thick 
Cdi_a;Mn2.Te layer. Dots with x < 0.035 are formed in 
Cdi_j,Znj,Te barriers with y ss 0.7. On the other hand, 
dots with X > 0.035 are formed in Zni_a;Mna;Te barri- 
ers with Mn concentration x kept nearly uniform in the 
barrier layers and in the dots. Such procedure ensures 
that Cdi-j^Mua^Te QDs with a high Mn content main- 
tain type-I band alignement as opposed to embedding 
the dots in ZnTe. Indeed, incorporation of a large den- 
sity of Mn ions into the Cdi_j.Mnj,Te dot leads to an 
increase of the QD band gap. Since the valence band 
offset in a CdTe /ZnTe heteroj unction is negligible, most 
of the hole confinement results from strain-induced band 
shifts. Increasing the Mn content decreases the lattice 
mismatch between Cdi_a;Mna;Te and ZnTe and thus di- 
minishes strain. As a result, above a certain value of x the 
dot potential for holes can become repulsive. Moreover 
this procedure excludes any effects related to interface 
interdiffusioiPSl and ensures a homogenous distribution 



of Mn ions in the dots. In the following, we refer to 
particular samples by their nominal Mn content x. 

We measure cw or time-resolved photoluminescence 
(PL) of single QDs excited with a 532 nm solid state 
laser or with a frequency doubled output beam of a op- 
tical parametric oscillator pumped with a Ti:Sapphire 
picosecond pulsed laser, respectively. Laser beam is fo- 
cused with a large numerical aperture microscope objec- 
tive. In order to access single dots, apertures with diam- 
eters of 500 nm are produced by spin-casting polystyrene 
beads and evaporating a 200 nm thick gold mask. The 
beads are subsequently lift-off by rinsing the sample in 
trichloroethylene and methanol. PL signal is collected 
with the same objective and analyzed with a CCD cam- 
era and a monochromator. Time-resolved detection is 
performed with a synchroscan streak camera with a tem- 
poral resolution of about 15 ps. For zero magnetic field 
measurements the sample is kept in a cold finger cryostat 
at 10 K. For the studies of the Zeeman effect, the sam- 
ple is immersed in pumped helium at 2 K in a split-coil 
cryostat providing fields up to 6 T in Faraday configura- 
tion. Two circular polarizations of the emitted light are 
recorded. 



IV. EXPERIMENTAL RESULTS 
A. Three relaxation stages 

In a previous reporPS, we demonstrated that the equi- 
librium between the Mn ions and the exciton (i.e., the 
EMP formation) is reached after the exciton has relaxed 
its spin and that establishing of the equilibrium during 
the exciton lifetime requires a concentration of Mn ions 
X > 0.1. Thus, for a sufficiently high x, the system 
is completely unrelaxed right after the photoexcitation, 
then the exciton adjusts its spin and then the EMP is 
formed. If the Mn density is not high enough, the exci- 
ton relaxes its spin, but the formation of the polaron is 
interrupted by recombination. In Section fll C| we showed 
that each of the three relaxation stages leads to a partic- 
ular shape of the PL spectrum. Now, we compare this 
prediction with the experiment. 

In Fig. [3] we show transient PL spectra at different 
time delays after the excitation pulse. For the sample 
with X = 0.01 (left panel), the spectrum is broad shortly 
after the excitation (top spectrum) and with increasing 
delay it narrows and becomes asymmetric without chang- 
ing much the central energy position (middle and bottom 
spectra) . The spectrum is noticeably steeper on the high 
than on the low energy side. This is a clear indication of 
the exciton spin relaxation (confront Fig. ITl. The com- 
plete equilibrium is not reached since the expected EMP 
form ation t ime for this Mn density is in the nanosecond 
rangcI^nSfilFor the sample with x = 0.2 (right panel), ex- 
citon spin relaxation is so efficient that the initial spec- 
tra recorded only 8 ps after photoexcitation (top panel) 
already exhibit a certain degree of asymmetry. It be- 
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FIG. 3. Transient PL spectra for a QD with x = 0.01 (inte- 
grated over 15 ps, left panel) and a; = 0.2 (integrated over 8 
ps, right panel) for three time delays At after the excitation 
pulse. 



comes even more clear as we look at larger delays (mid- 
dle panel). Simultaneously, the spectrum undergoes a 
redshift reflecting lowering of the exciton energy upon 
EMP formation, which for this QD develops in t/ = 100 
ps. After reaching the equilibrium, the Mn ions are po- 
larized giving rise to a magnetization resulting from the 
effective field B^x imposed by the exciton. Consequently, 
the spectrum becomes symmetric and Gaussian with its 
FWHM reflecting the magnetization fluctuations within 
the polaron. Here, we do not attempt fltting of the PL 
line shapes with Eqs. (14 1 and (19 1, since these spectra 



do not pertain to any steady states, but rather reflect 
mixed stages of a dynamical process. However, it is clear 
that the distinction of speciflc relaxation stages is correct. 

Another qualitative agreement between our line shape 
theory and the experiment is demonstrated in Fig. [4] 
where we compare measured and calculated PL spectra 
at various temperatures for two dots: with x — 0.035 
(left panel) and with x = 0.2 (right panel). In the for- 
mer case, we observe two transitions, one related to the 
neutral exciton (X°) and the other to a charged exci- 
ton {X*), both of them asymmetric as a consequence of 
the spin-relaxed exciton population. For this Mn con- 
centration, the exciton spin relaxation time is on the or- 
der of 10 ps^ZI — much shorter than the recombination 
time, which amounts to about 300 ps.^S We thus observe 
that these spectra reflect a steady state with a relaxed 
exciton spin. As the temperature is increased, the im- 
balance between occupation probabilities of the spin-up 
and spin-down excitons states becomes less pronounced 
and the PL spectrum evolves into a symmetric Gaussian 
(see Eq. 19 1. Simultaneously, the transitions red-shift as 
a result of the shrinkage of the band gap. On the other 



hand, the spectrum of the dot with x = 0.2 at low tem- 
perature is approximately a Gaussian as a result of the 
total equilibrium attained here during the exciton life- 
time (see Fig. Inland the discussion above). Note that in 
these time-integrated spectra an additional broadening 
due to the transient energy shift reflecting the EMP for- 
mation can be present. As the temperature is increased, 
the magnetization of the EMP is diminished and conse- 
quently the exciton recombination energy is blueshifted. 
However, the excitons remain spin-relaxed and thus with 
increasing the temperature the PL spectrum acquires an 
asymmetry. As the temperature is further increased the 
spectrum broadens, redshifts, and becomes symmetric as 
for the X ~ 0.035 sample. 
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FIG. 4. PL spectra measured for a QD with x — 0.035 (a) and 
X = 0.2 (b) at various temperatures compared with calculated 
ones in, respectively, (c) and (d). For the calculation in (c), 
spin relaxed excitons are assumed. For the calculation in (d), 
we assume the EMP formation resulting from the effective 
field Bex- The EMP is suppressed with increasing tempera- 
ture. The spectra are calculated for the same temperatures 
as in the experiment and coded with the same colors. 



A direct comparison between our transition line model 
and the experiment is shown in Fig. p[a), where we 
present PL spectra for a QD with x = 0.035 immersed 
in liquid helium at 2 K, in magnetic flelds up to 6 T, 
positive (negative) sign denoting ct"*" (a^ ) polarization of 
detection. The inset shows the close-up of the spectra for 
the smallest flelds, clearly demonstrating the line shape 
anisotropy. The asymmetry of the spectrum is more pro- 
nounced than for the spectrum for a dot with the same 
Mn content shown in Fig. [4] where lowest bath tempera- 
ture was 9 K. Above 0.5 T, the signal in a^ polarization 
disappears reflecting the vanishing of the occupation fac- 
tors introduced in Eq. (19 1. Points and lines denote the 



measured and fltted spectra, respectively. In the fltting, 
we assume completely spin relaxed excitons (see expres- 
sion ( |19| ). The fltting is performed simultaneously to the 
whole set of 27 spectra from -0.5 to +6 T with a step of 
0.25 T (for the sake of visibility, in Fig. [Elwe show every 
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2nd spectrum). We keep constant the model parameters 
unaffected by the magnetic field, i.e., A^Mn, -^o^^j and 
Tx, while letting the transition amplitude and Mn spin 
temperature evolve with the magnetic field. 
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concentration XeS = 0.024 and thus, in the studied dot 
we expect to have A^o^ • Xes = 96 paramagnetic Mn ions. 
We thus find the fitted value in excellent agreement with 
this expectation, since the number fluctuations are on 
the order of \/NMn ~ 10. 

In Fig. IsFb), we show the fitted, magnetic field depen- 
dent Mn spin temperatures. For magnetic fields smaller 
than 0.5 T, the Mn spin system is efficiently heated by 
the energy transfer from the hot, photocr eated carriers 
by simultaneous carrier-ion spin flip-flops.H3Sl In that 
case Tx controls the PL lineshape and the Mn spin tem- 
perature becomes irrelevant as long as its value is higher 
than about 10 K. As the magnetic field is increased, a 
cooling process occurs as the flip-flops become energeti- 
cally forbidden when the exchange driven carrier Zeeman 
splitting becomes larger than Mn spin splitting governed 
by the Mn g- factor of 2.13ZES! As a consequence, the Mn 
spin temperatures decrease from about 15 down to 8 K. 
For magnetic flelds above about 3 T, the spin tempera- 
ture is stabilized as the Mn ions attain a thermal equilib- 
rium with the lattice. We thus conclude that the lattice 
temperature is about 8 K. It remains a few kelvin larger 
than the bath temperature of 2 K. We expect that this 
elevated temperature originates from diffusive heat trans- 
port via the LO phonons cascaded down in the energy re- 
laxation process of nonresonantly excited carriers.*^ The 
fitted value of Tx = 9.9K is slightly higher than the 
lattice temperature, which is a consequence of the time- 
integrated nature of this PL measurement (see discussion 
in Sec. II C). In the inset to Fig. Islb), we plot the total 
transition intensities as a function of the magnetic field 
obtained by integrating the fitted spectra. We observe 
an increase by about a factor of 10 between and 4 T. 
We attribute this effect to an increased PL yield resulting 
from blocking of the excita tion tran sfer from the QDs to 
the internal Mn transitions.'SSESIIIl 



FIG. 5. (a) Measured (points) and fitted (lines) PL spectra for 
a QD with x = 0.035. The fitting is performed for the entire 
set of spectra simultaneously yielding parameters A^ivin = 99, 
A^o^ = 3200. The inset shows a close-up of the small field 
spectra, (b) Magnetic field dependent Mn spin temperature 
obtained from the fitting. Inset: the fitted, magnetic field 
dependent integrated intensity. 

The accuracy of the fit is rather excellent. Assuming a 
lens shape QD, the obtained localization volume NqV ^ 
3985 corresponds to a dot with 19 nm in diameter and 2 
nm in height, which is a typical size.l^The fitted number 
of paramagnetic Mn ions A^Mn = 99 should be analyzed 
carefully, since this model parameter not only includes 
those ions which do not have a Mn as a nearest neighbor, 
but also a combined effect of further neighbors leading 
to a suppressed paramagnetic response (see discussion in 
Section II A I . Moreover, statistical fluctuations of the Mn 
number have to be considered. Following the discussion 
presented in Section [IIA| we can calculate the number of 
Mn ions without a Mn as a nearest neighbor. For dots 
with X = 0.035 grown on Cdi_j,Znj,Te the effective Mn 



B. Transition energy and line\vidth 

In this Section, we analyze the transition energy (en- 
ergy of the PL maximum) and the transition linewidth 
(FWHM) as a function of the magnetic fleld for 
Cdi_^Mn^Te QDs with x ranging from 0.01 to 0.2. The 
objective of this analysis is to demonstrate how the spec- 
troscopic features change with the Mn content, to retrieve 
the morphological parameters from the PL spectra, and 
to discuss what parameters control the particular features 
at different Mn densities. 

In Fig. p] we show excitonic Zeeman shifts (top) and 
linewidths (bottom) for three Cdi-ajMn^jTe dots: with 
nominal x = 0.01 (Fig. [6la,b)), with a; = 0.05 (Fig. 
l6|c,d)), and with x = 0.1 (Fig. [6Fe,f)). Measured Zee- 
man shifts are fitted with Brillouin functions according 
to Eq. ((lel). The Hnewidths are fitted with Eq. p7|. 



For a given dot, the fitting is performed simultaneously 
for the two sets of data. Analogous fits are obtained for 
the two remaining samples, with x = 0.035 and x = 0.02. 
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FIG. 6. Zeeman shifts (top) and linewidths (bottom) for exciton transitions from Cdi_a;Mna;Te QDs with intended Mn contents 
of X = 0.01 (a,b) X = 0.05 (c,d), and x = 0.1 (e,f). Points are experimental values, lines are Zeeman shifts and linewidths fitted 



according to Eqs. (16 1 and (171, respectively. Fitting is performed simultaneously for the energy shift and linewidth for each 



particular dot. For the dots with x = 0.01 and with x = 0.1, the data points for low magnetic fields are skipped from the fitting 
procedure (see dashed lines and the discussion in the text) 



The fitting parameters are A^Mn, -^o^j ^^nd Mn spin tem- 
perature T. This temperature is assumed here as field 
independent and thus the fitting range for samples with 
X < 0.05 are adjusted by omitting the data points below 
2 T, where heating of the Mn spins occurs, as discussed 
above. Moreover, an anomalous non-monotonicity of the 
transition linewidth is observed for dots with smallest Mn 
concentrations (see Fig. l6jb)) which cannothe accounted 
for by the heating effects. For samples with x > 0.05 we 
do not expect any significant variations of the Mn spin 
temperature with the magnetic field, since for such large 
Mn densities the heat transferred from the hot photo- 
carriers is efficiently dissipated via the spin-lattice relax- 
ation, which is e nhanc ed in the presence of stronger spin- 
spin interactions.l^fiEl] However, for the samples in which 
the complete EMP formation occurs during the exciton 
lifetime {x = 0.1 and x = 0.2), the transition redshift sig- 
nificantly affects the time-integrated PL linewidth (see 
data points below 1 T in Fig. l6jf)), and hence data 
points below 1 T are omitted. For these samples, we fit 
the functions from Eqs. ( 16 1 and (17) assuming that the 
Mn ions experience the total field B — Bo + B^x, where 
Bq is the external magnetic field and B^^ is given by Eq. 



(201 



Fit accuracy is very good. Obtained fitting parameters 



for the samples identified by x are collected in Table 1. It 
shows that the number of paramagnetic ions in a dot can 
be as high as about 200 provided that the dot volume is 
large enough. Above x = 0.05, the Nmh decreases. This 
is expected since the concentration of Mn cations with- 
out a Mn as a nearest neighboralso decreases. Obtained 
values of the exciton localization volume are expressed in 
the number of cation sites and vary from about 1000 to 
5500. These numbers agree very well with the QD sizes 
obtained in the atomic force microscopy (AFM) study of 
an uncapped layer of Cdi-a^Mn^^Te dots with x — 0.03 
grown on Cdi_j,ZnyTe with y « 0.8. For lens shaped 
QDs with a height of 2 nm, our range of volumes corre- 
spond to diameters from about 10 to 22 nm. For the dots 
with the largest Mn density of a: = 0.2, the fitted localiza- 
tion volume is particularly small, which could be due to 
exciton autolocalization by the spatially inhomogenous 
EMP potential.^' '^^ Obtained Mn spin temperatures are 
again a few kelvin higher than the bath temperature and 
as above, we conclude that this originates from heating of 
the lattice itself by the nonresonant excitation. In Table 
1, we also compare the obtained values of Nmu with those 
calculated N^^ = NoVxeS, where A^o^ is the fitted local- 
ization volumes and XcS = xSo{x)/S, where Sq{x) is the 
effective spin saturation value applied to parameterize 
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the Brillouin functiorP^BSl xhe agreement is very good 
for smallest Mn concentrations and in all cases the cal- 
culated values lie within three standard deviations from 
the fitted ones. 



X 


A''m„ 


NoV 


T 


N^n 


0.01 


24 


2970 


3.1 K 


26 


0.035 


41 


2040 


4.2 K 


46 


0.05 


183 


5270 


8.0 K 


150 


0.1 


69 


2230 


7.0 K 


89 


0.2 


62 


1120 


5.4 K 


48 



TABLE I. The parameters obtained by simultaneous fitting 
of the transition energy and linewidth for five QD with differ- 
ent nominal Mn compositions. Exciton localization volume 
NqV is expressed in cation sites. Fitted numbers of paramag- 
netic Mn ions are compared with the ones calculated (A'^Mn) 
assuming a suppression of paramagnetic response by the near- 
est neighbors. 

In Fig. IrFa), we demonstrate Zeeman shifts (points) 
measured at 5 T with the sample at bath temperature of 
2 K averaged over 3-10 dots for each Mn concentration x. 
These values are compared (brown broken line) to Zee- 
man shifts calculated using Eq. ( 16 ) for T = 7 K and 



replacing the ratio Nmii/NqV by XeS calculated above. 
Analogously to fitting performed in Fig. |6] for a; > 0.05 
in the calculation we put _B = 5 T -I- Bex accounting for 
the EMP formation and B^x is calculated for an average 
exciton confinement volume spanning over 2000 cation 
sites. Taking into account that the experimental values 
correspond to a rather small statistical ensemble, the av- 
erage values put into the Zeeman shifts calculation re- 
produce the overall trend very well. 
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FIG. 7. (a) Exciton Zeeman shifts and (b) transition 
linewidths averaged over a set of dots with different Mn con- 
tent X. Broken line in (a) shows calculated shifts assuming a 
statistical distribution of Mn ions (see text). 

Fig. m[^b) shows transition linewidths averaged over 
sets of 10-20 dots for each of the samples with different 



Mn contents, measured at bath temperature of 10 K and 
zero external magnetic field. The linewidth monotoni- 
cally increases with increasing Mn concentration. How- 
ever, this increase is driven by different mechanisms, de- 
pending on X. For dots with x < 0.035, narrow linewidths 
result from the narrowing due to exciton spin relaxation 
despite increased Mn spin temperature. On the other 
hand, for dots with x > 0.035 the transitions are broad- 
ened due to the transient shifts reflecting the EMP for- 
mation (still influencing the a; < 0.1 dots even if the 
full equilibrium is not attained during the exciton life- 
time). Hence, the increase of transition linewidth cannot 
be directly linked to the enhanced magnetization fluctua- 
tions, since at different Mn concentrations various mech- 
anisms influence its value. In fact, Eq. [iTJshows that the 
linewidth resulting from the magnetization fluctuations 
scales as ^/NMn/iNoV)'^ and thus for dots at the same 
stage of the equilibration in the exciton-Mn spin system, 
we would expect the linewidth to strongly decrease for 
dots with highest Mn contents. 



V. DISCUSSION 

Experimental results presented in Figs. |3] |4l and [5] 
unambiguously demonstrate spin relaxation of excitons 
in the studied Cdi-^^Mn^Te QDs. This is an important 
conclusion, since a vast majority of the studies of spin 
relaxation in QDs focuses on spin relaxation of single 
carriers. A natural question arises about the mechanism 
of the exciton spin relaxation. We propose two types 
of mechanisms, differing in the role played by the sp-d 
exchange interaction. 

In the flrst mechanism, the role of the Mn ions is static 
and the role of the thermal reservoir is played by phonons. 
The exchange interaction provides a splitting of the |±1) 
exciton states. Momentary fluctuations of the QD mag- 
netization results in exchange splittings, which in turn in- 
fluence the transition linewidth of a time-integrated spec- 
trum. The exchange splitting can be approximated from 
the transition linewidth taking into account the line nar- 
rowing due to exciton spin relaxation or broadening due 
to transient energy shift during EMP formation. E.g., 
from the PL spectra for a dot with x = 0.035 presented 
in Fig. [5] we can estimate the exchange induced split- 
ting of the exciton state to be equivalent to a direct Zee- 
man splitting of an exciton in CdTe dot in a magnetic 
fleld of about 15 T. In order for the spin-flip transition 
due to phonon scattering to occur, the Zeeman split sub- 
levels need to have a mixed spin character. -=^^ For non- 
magnetic QDs, the intermixing can be provided by e.g., 
the combined effect of the strain-driven heavy-light hole 
mixing and the electron-hole exchange interaction. The 
mixing allows a phonon-assisted spin-flip transition and 
the resulting relaxation rates between the Zeeman split 
bright exciton states scale as^-: l/Ti ~ s^{Eih/do)'^B^, 
where s is the sound velocity and Eh^/Sq is the size- 
dependent ratio of the heavy-light hole splitting and 
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the isotropic exchange sphtting (see Section II B ) . For 



InAs/GaAs QDs at 15 T, this mechanism yields spin re- 
laxation times of the order of a few nanoseconds, larger 
for smaller dots. For a CdTe dot, these times can be 
further decreased as a result of the sound velocity be- 
ing smaller by about 40% allowing to obtain values com- 
parable to the exciton lifetime and in the same order 
of magnitude as measured experimentally.'^^ Moreover, 
the relaxation between the bright and dark states is ex- 
pected to be even more efficient.'^ Another mixing mech- 
anism allowing for the phonon-assisted spin-flip is due to 
the spin-orbit coupling. This mechanism combined with 
piezo-electric spin-phonon coupling was sh own to deter- 
mine the relaxation rates of single electrongSEH confined 
in QDs and to scale as l/Ti r^ B^ , providing an even 
more efficient process. Remarkably, in Cdi-ajMn^jTe QDs 
the s-d exchange interaction provides another important 
mechanism for mixing of the exciton states, namely via 
electron-Mn ion spin flip-flops (see Hamiltonian ^ and 
the discussion in Section II B I . Therefore, due to the com- 



bined action of the carrier-Mn exchange, electron-hole ex- 
change, and the heavy-light hole mixing, the spin-up and 
spin-down exciton states are substantially mixed which 
in itself allows to expect short relaxation times. It also 
makes a rigid theoretical modeling challenging. 

The second spin relaxation mechanism involves a dy- 
namical role of the Mn ions. The sp-d exchange is no 
longer just a source of an effective quasi-static magnetic 
field giving rise to a splitting of the Kramers doublet, 
but may act as a scattering mechanism due to the dy- 
namics of Mn spins caused by their short-range inter- 
actions (both isotropic and anisotropic superexchange) . 
Such a role was invoked for spin relaxation of electrons in 
CdTe/Cdi_^Mn:rTe quantum wells.'^Sl Proper description 
of the carrier (electron, hole, and exciton) spin relaxation 
in semimagnetic QDs obviously requires more theoretical 
investigations. 

The model of the PL spectrum presented here is able 
to reproduce experimental data in a relatively wide range 
of external parameters such as magnetic field and tem- 
perature. The main effect not accounted for is the co- 
existence of different charge states in the single dot PL 
spectrum. As shown for CdTe QDs, even at small ex- 
citation powers all the s-shell transitions appear simul- 
taneously in the spectrum, with a universal sequence 
E{X^) > E{X+) > E{X-) > E{2X), where Ej x) is 
the transition energy of the excitonic complex t^BSESI 
For Cdi_j;Mnj;Te dots with a sufficiently low Mn con- 
tent it is possible to identify the neutral exciton and a 
charged exciton (see e.g.. Fig. El, but due to the mag- 
netization fluctuations, the broadening of the transitions 
result in charge exciton recombinations merged together 
with the biexciton transition. For dots with Mn con- 
centration above 0.1, the transition linewidth ultimately 
precludes identification of the charge states. Our the- 
oretical considerations assume a neutral QD and thus 
the PL spectrum model and discussion of the magnetic 
field dependence of the Zeeman shifts and linewidths con- 



cern exclusively the recombination of a neutral exciton. 
However, emission from a charged dot can exhibit a very 
different behavior in magnetic field, especially with re- 
spect to the development of carrier-Mn ion equilibrium. 
A positively charged dot, even with a very small Mn con- 
centration, can reach the equilibrium due to the presence 
of the additional hole not taking part in the recombina- 
tion and dwelling in the QD for a time much longer than 
both the exciton lifetime and EMP formation. This hole 
may develop a net magnetization due to the hole-Mn ion 
exchange interaction. For a dominant X^ recombination 
we therefore expect the signatures of equilibrium such as 
the transition blueshift with increasing temperature or 
decreased Zeeman shift to appear even when the Mn den- 
sity is low and EMP formation time is large. For a high 
Mn content, the picture becomes more complicated since 
both the single carrier and the exciton may develop the 
magnetization. However, while photocreation of an X'^ 
turns on the formation process, the excitation of the X+ 
is rather expected to at least diminish the magnetization, 
since the two holes form a singlet states. Nevertheless, 
we expect to see different Zeeman shifts for a neutral and 
charged dot. 



VI. CONCLUSIONS 

We have investigated the properties of the PL spec- 
trum of semimagnetic, Cdi_a;Mna;Te QDs. We found 
that in the process of reaching a full thermal equilibrium 
in the exciton-Mn ion spin system, the PL spectrum un- 
dergoes changes in energy and linewidth. In particular, 
the equilibration process involves an intermediate stage, 
in which the exciton spin is relaxed, while the full exciton- 
Mn ion system is not yet in equilibrium. This particu- 
lar stage gives rise to a specific shape of the spectrum, 
with a substantially narrowed and asymmetric transition 
line. We have developed a theoretical model, allowing to 
analytically compute the PL spectra for the unrelaxed, 
exciton spin-relaxed, and fully relaxed system. We found 
an excellent agreement between the model calculations 
and the experiment, and the comparison allowed to ac- 
cess such parameters as the exciton localization volume, 
the number of (approximately) paramagnetic Mn ions, 
and their temperature. The volumes evaluated from the 
experiment remain in agreement with AFM studies. The 
numbers of paramagnetic Mn ions in dots with different 
Mn densities agree well with those calculated assuming 
statistical distribution of the impurities in the cation sub- 
lattice. Furthermore, we have found that the heating of 
the Mn ions by the photocreated carriers is not effective 
for relatively large Mn densities although the Mn ions are 
still hotter than the bath. We have presented a compre- 
hensive study of the PL energy and transition linewidth 
for Cdi-^cMuj^Te quantum dots in a broad range of Mn 
concentrations and discussed the processes controlling 
these parameters for various Mn densities. 

This research was supported by Polish Min- 



12 



istry of Science and Higher Education grant no. 
2011/01/B/ST3/022287 and by European Union within 
European Regional Development Fund, through Innova- 
tive Economy grant (POIG.OI. 01. 02-00-008/08). L.C. 
acknowledges support from the "FunDMS" Advanced 
Grant of the ERG within the "Ideas" 7th Framework 
Programme. We thank Tomasz Dietl and Jacek Kossut 
for discussions. 



Appendix A: Effective volume determined by 
interaction of the exciton with the Mn spins 

Let us define the ratio of hole and electron volume as: 



i? = 



Vh 



(Al) 



and assume that 1/2 < i?< 2. We will now consider two 
cases: (i) Mn spins located only within the QD volume 
(which we identify here with the Vg) and (ii) Mn spins 
uniformly distributed in the QD and the surrounding bar- 
rier. Taking into account that in Gdi-^^Mn^^Te |/3/a| =4, 
the effective volume defined in Eq. ( 12 1 is then given by: 



V = 



hVh 



R 



(A2) 



in case (i) and V — Vh in case (ii). Hence, for the assumed 
range of values for R, we get that in the former case 
l()/^-Vh>V>b/&-Vh^Vh. 

Note that there is an alternative definition of the ef- 
fective volume, in which we take the Hamiltonian from 



Eq. (11) and use it to derive the expression for cje anal- 



ogous to the one from Eq. (15 1: 



{a-PY 2 



4K.2 



-o-s 



(A3) 



defining in this way the volume ¥„. A short calculation 
shows that in case (i) we have Vo- — V . while in case (ii) 
we have: 



K = 



5 Vh 



4./^T^ ' 

2 ' 16 



(A4) 



which gives 1.01 ■Va>V> 0.98 • Va- ~ Vh- This shows 
that both definitions of the effective volume are, with 
reasonable accuracy, consistent with each other and with 
the statement that the effective volume obtained from 
the fitting (Figs. (5|6)) corresponds to the confinement 
volume of the hole wavefunction. 



* Corresponding author: lukasz.klopotowski@ifpan.edu.pl 
^ P. Michler, ed., Single Semiconductor Quantum Dots 

(Springer- Verlag, 2009). 
^ Y. Arakawa and H. Sakaki, Appl. Phys. Lett., 40, 939 

(1982). 
^ P. Michler, A. Kiraz, C. Becher, W. V. Schoenfeld, P. M. 

Petroff, L. Zhang, E. Hu, and A. Imamoglu, Science, 290, 

2292 (2000). 



* A. Barenco, D. Deutsch, A. Ekert, and R. Jozsa, Phys. 
Rev. Lett., 74, 4083 (1995) 



* D. Loss and D. P. DiVincenzo, Phys. Rev. A, 57, 120 
(1998). 

® M. Kroutvar, Y. Duconiniun, D. Heiss, M. Bichler, 

D. Schuh, G. Abstreiter, and J. J. Finley, Nature, 432, 81 

(2004). 
'^ R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 

and L. M. K. Vandersypen, Rev. Mod. Phys., 79, 1217 

(2007). 

* R.-B. Liu, W. Yao, and L. J. Sham, Adv. Phys., 59, 703 
(2010). 

® J. Houel, A. V. Kuhlmann, L. Greuter, F. Xue, M. Poggio, 
B. D. Gerardot, P. A. Dalgarno, A. Badolato, P. M. Petroff, 
A. Ludwig, D. Reuter, A. D. Wieck, and R. J. Warburton, 



Phys. Rev. Lett., 108, 107401 (2012) 



J. Fernandez-Rossier and L. Brey, Phys. Rev. Lett., 93, 

117201 (2004). 
" A. O. Govorov, Phys. Rev. B, 72, 075359 (2005). 
^^ R. M. Abolfath, P. Hawrylak, and L Zutic, Phys. Rev. 

Lett., 98, 207203 (2007). 
^^ S. Mackowski, G. Prechtl, W. Heiss, F. V. Kyrychenko, 

G. Karczewski, and J. Kossut, Phys. Rev. B, 69, 205325 



(2004). 
^^ T. Gurung, S. Mackowski, G. Karczewski, H. E. Jackson, 



and L. M. Smith, |App. Phys. Le tt., 93, 15 3114 (20 0JT[ 



Y. Lyanda-Geller, T. L. Reinecke, and G. Bacher, [Journal 



t 



of A pplied Ph ysics, 111, 09 3705 (2012) 
^^ A. Kudelski, A. Lemaitre, A. Miard, P. Voisin, T. C. M. 
Graham, R. J. Warburton, and O. Krebs, Phys. Rev. Lett"^ 
99,247209(2007), 

J. A. Gaj and J. Kossut, eds.. Introduction to the Physics of 
Diluted Magnetic Semiconductors (Springer- Verlag Berlin 
Heidelberg, 2010). 

A. A. Maksimov, G. Bacher, A. McDonald, V. D. Ku- 
lakovskii, A. Forchel, C. R. Becker, G. Landwehr, and 
L. W. Molenkamp, Phys. Rev. B, 62, R7767 (2000) 
G. Bacher, A. A. Maksimov, H. Schomig, V. D. Ku- 
lakovskii, M. K. Welsch, A. Forchel, P. S. Dorozhkin, A. V. 
Chernenko, S. Lee, M. Dobrowolska, and J. K. Furdyna, 
Phys. Rev. Lett., 89, 127201 (2002) 

A. Hundt, J. Puis, and F. Henneberger, Phys. Rev. B, 69, 
121309 (2004) 



L. Besombes, K. Kheng, L. Marsal, and H. Mariette, Phys. 
Rev B, 63, 155307 (2001)| 



^^ D. Gammon, E. S. Snow, B. V. Shanabrook, D. S. Katzer, 
and D. Park, Science, 273, 87 (1996). 

^^ P. Wojnar, J. Suffczynski, K. Kowalik, A. Golnik, G. Kar- 
czewski, and J. Kossut, Phys. Rev. B, 75, 155301 (2007)| 

^'^ J. Seufert, G. Bacher, M. Scheibner, A. Forchel, S. Lee, 
M. Dobrowolska, and J. K. Furdyna, |Phys. Rev. LettT^ 
88, 02740 2 (2001) , 

^^ R. Beaulac, L. Schneider, P. I. Archer, G. Bacher, and 
D. R. Gamelin, Science, 325, 973 (2009). 



13 



I. R. Sellers, R. Oszwaldowski, V. R. Whiteside, M. Egin- 
ligil, A. Petrou, I. Zutic, W.-C. Chou, W. C. Fan, A. G. 
Petukh ov, S. J. Kim, A. N . Cartwright, and B. D. Mc- 
Combe, Phys. Rev. B, 82, 195320 (2010) 
L. Klopotowski, L. Cywinski, P. Wojnar, V. Volio- 
tis, K. Fronc, T. Kazimierczuk, A. Golnik, M. Ravaro, 
R. Grousson, G. Karczewski, and T. Wojtowicz, Phys. 



Rev. B, 83, 081306 (2011) 



^° J. M. Pientka, R. Osz waldowski, A. G. Pet ukhov, J. E. 

Han, and I. Zutic, Ph ys. Rev. B, 86, 161403 (2012) 
^^ K. V. Kavokin, I. A. Merkulov, D. RTYakovlev, W. Ossau, 

and G. Landwehr, Phys. Rev. B, 60, 16499 (1999) 
^° T. Dietl, P. Peyla, W. Grieshaber, and Y. Merle 

d'Aubigne, Phys. Rev. Lett., 74, 474 (1995). 
31 D. R. Yakovlev, Acta Phys. Pol. A, 90, 703 (1996). 
3^ J. A. Gaj, R. Planel, and G. Fishman, Solid State Com- 

mun., 29, 435 (1979). 
33 T. Dietl and J. Spalek, Phys. Rev. B, 28, 1548 (1983). 
3* R. Brazis and J. Kossut, Solid State Coniniunications, 122, 

73 (2002). 
3^ L. Besombes, Y. Leger, L. Maingault, D. Ferrand, H. Ma- 

riette, and J. Cibert, [Pivs. Rev. Lett., 93, 207403 (2004)| 
3® Y. Leger, L. Besombes, L. Maingault, and H. Mariette, 

Phys. Rev. B, 76, 045331 (2007) 
3^" M. Goryca, P. Plochocka, T. Kazimierczuk, P. Wojnar, 

G. Karczew ski, J. A. Gaj, M. Potemski, and P. Kossacki, 

Phys. Rev. B , 82, 165323 (2010||_ 



L. Cywihski, Phys. Rev. B, 82 , 075321 (2010) 1 
T. Smolenski, T. Kazimierczuk, M. Goryca, T. Jakubczyk, 
Lukasz Klopotowski, Lukasz Cywinski, P. Wojnar, A. Gol- 
nik, and P. Kossacki, Physical Review B, 86, 241305 
(2012). 

M. Bayer, G. Ortner, O. Stern, A. Kuther, A. A. Gor- 
bunov, A. Forchel, P. Hawrylak, S. Fafard, K. Hinzer, T. L. 
Reinecke, S. N. Walck, J. P. Reithmaier, F. Klopf, and 
F. Schafer, [PFys. Rev. B, 65, 195315 (2002) 
P. Wojnar, J. Suffczynski, K. Kowalik, A. Golnik, 



M. Aleszkiewicz, G. Karczewski, and J. Kossut, Nanotech- 

nology, 19, 235403 (2008). 
'^^ L. Besombes, Y. Leger, L. Maingault, and H. Mariette, 

Journal of Applied Physics, 101, 081713 (2007) 
■13 F. Tinjod, B. Gilles, S. Moehl, K. Kheng, and H. Mariette, 

Appl. Phys. Lett., 82, 4340 (2003) 
"i* P. Wojnar, C. Bougerol, E. Bellet-Amalric, L. Besombes, 

H. Mariette, and H. Boukari, J. Cryst. Growth, 335, 28 

(2010). 
^^ W. Grieshaber, A. Haury, J. Cibert, Y. M. d'Aubigne, 

A. Wasiela, and J. A. Gaj, Phys. Rev. B, 53, 4891 (1996). 
^^ L. Klopotowski, V. Voliotis, A. Kudelski, A. L Tar- 

takovskii, P. Wojnar, K. Fronc, R. Grousson, O. Krebs, 

M. S. Skolnick, G. Karczewski, and T. Wojtowicz, Phys. 

Rev. B, 83, 155319 (2011). 
■'■^ B. Konig, L A. Merkulov, D. R. Yakovlev, W. Ossau, S. M. 

Ryabchenko, M. Kutrowski, T. Wojtowicz, G. Karczewski, 

and J. Kossut, Phys. Rev. B, 61, 16870 (2000). 
*^ T. Clement, D. Ferrand, L. Besombes, H. Boukari, and 

H. Mariette, Phys. Rev. B, 81, 155328 (2010). 
^^ A. Hundt J. Puis, A. V. Akimov, Y. H. Fan , and F. Hen- 

neberger, |Phys. Rev. B, 72, 033304 (2005)| 
^° C. Kim, M. Kim, S. Lee, J. Kossut, J. Furdyna, and 

M. Dobrowolska, J. Cryst. Growth, 214/215, 395 (2000). 
^1 M. Nawrocki, Y. G. Rubo, J. P. Lascaray, and D. Coquil- 

lat, Phys. Rev. B, 52, R2241 (1995). 

A. V. Khaetskii and Y. V. Nazarov, Physical Review B, 

64, 125316 (2001). 



52 



53 E. Tsitsishvili, R. v. Baltz, and H. Kalt, |Phys. Rev. B, 67, 
205330 (2003) 

="* K. Roszak, V. M. Axt, T. Kulm, and P. Machnikowski, 
Phys. Rev. B, 76, 195324 (2007) 

55 R. Akimoto, K. Ando, F. Sasaki, S. Kobayashi, and 
T. Tani, Phys. Rev. B, 56, 9726 (1997). 

5® T. Kazimierczuk, T. Smolenski, M. Goryca, L. Klopo- 
towski, P. Wojnar, K. Fronc, A. Golnik, M. Nawrocki, J. A. 
Gaj, and P. Kossacki, Phys . Rev. B, 84, 165319 (2011)| 



